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Abstract. The dissolution time (tdis) of clusters in a tidal field does not scale with the "classical" 
expression for the relaxation time. First, the scaling with N , and hence cluster mass, is shallower 
due to the finite escape time of stars. Secondly, the cluster half-mass radius is of little importance. 
This is due to a balance between the relative tidal field strength and internal relaxation, which 
have an opposite effect on tdis, but of similar magnitude. When external perturbations, such 
as encounters with giant molecular clouds (GMC) are important, tdis for an individual cluster 
depends strongly on radius. The mean dissolution time for a population of clusters, however, 
scales in the same way with mass as for the tidal field, due to the weak dependence of radius 
on mass. The environmental parameters that determine tdis are the tidal field strength and the 
density of molecular gas. We compare the empirically derived tdis of clusters in six galaxies 
to theoretical predictions and argue that encounters with GMCs are the dominant destruction 
mechanism. Finally, we discuss a number of pitfalls in the derivations of tdis from observations, 
such as incompleteness, with the cluster system of the SMC as particular example. 

Keywords, globular clusters: general, open clusters and associations: general, stellar dynamics, 
methods: n-body simulations 



1. Theoretical predictions of cluster dissolution 

1.1. Dynamical evolution in a tidal field 
Simulations of star clusters dissolving in a tidal field have shown that the dissolu- 



tion t ime jt ^At) scales with the relaxation time (trei) as tj^'^ oc t^^y^ (jBaumgardtl 12001 



Baumgardt fc Makinoll2003l ) . This non-line ar dependence on trei is due to the finite escape 



time through one of the Lagrange points ( Fukushige fc Heggie 2000() . The dependence 
on or cluster m ass (Mc), can be approximated as t*,'^ o c Mc°'^^, which is accurate 
for 10^ < TV < 10^ (|Lamers. Gieles fc Portegies ZwartI [20051) . The half-mass radius (rh) 
of the cluster does not enter in the results, since it is assumed that clusters are initially 
"Roche lobe" filling, which implies Mc oc rj^, i.e. a constant crossing time. 

The assumption of Roche lobe filling clusters is computationally attractive since it 
avoids having rh as an extra parameter. However, observations of (young) extra-galactic 
star clusters show that the dependence of rh on Mc and galactocentric distance {Rq) 
is considerably weaker (rh oc M^-^ Rq^) than the Roche lobe filling relation (rh oc 
M^/^Rq^^) ()Larsenl 12004 IScheeomaker. Haas. Gieles et a/.ll2007h . implying that mas- 
si ve clusters at large Rg are i nitially underfilling their Roche lobe. 

Gieles fc Baumgardt! ( 2007 ) simulated clusters with varying initial rh in a tidal field to 
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Figure 1. Half-mass time as found from A^-body simulations of clusters dissolving in a tidal 
field. The filled circles represent clusters that initially fill their Roche lobe. The open squares 
are the results of runs where rh was seven times smaller. 



quantify the importance of rh . Figure [T] shows the results of t for two sets of clusters 
with different initial rh. The filled circles are for clusters that started tidally limited and 
the open squares are for runs where the initial rh was a factor seven smaller. The difference 
in tdis are within a factor two, while the "classical" expression of t-^ci predicts a factor 
73/2 ^ 20. The reason that idis depends so little on rh can be understood intuitively: for 
smaller clusters the tidal field is less important, but the dynamical evolution is faster. 
These effects happen to balance and result in almost no dependence on rh. The crossing of 
the lines around N ~ 10^ implies that for globular clusters is completely independent 
ofri,. 

This somewhat surprising result means th at we can use the rh independent results for 
tdis of tidally limited clusters (Baumgardt fc Making 2003f l as a general result for f^'^ for 
clusters of different rh: 
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From this it follows that a cluster with = IC* Mq in the solar neighbourhood would 
dissolve in approximately 8 Gyr due to tidal field. This is much l onger than the empirically 



derived value of 1.3 Gyr (jLamers. Gieles, Bastian. et aLl 120051 ). implying that there are 



additional disruptive effects that shorten the life-time of clusters. 



1.2. External perturbations: disruption by giant molecular clouds 

It has long been suspected t hat encounters with giant inolecular clouds (GMCs) shorten 

the lif e-times of clusters (e.g. van den Bergh fc McClurellQSd ). Gieles. Portegies Zwart. Baumgardt. et al. 



studied this effect using A''-body simulations and found that tdis due to GMC en- 
counters (^di^^) can be expressed in cluster properties and average molecular gas density 
(/On) as 



iSf ^ _ „ . / 0.03A/epc-n / M, \ / 3.75pc 
Gyr pn )Wmq) I rh ' " ^'''^ 
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Table 1. Columns 1-3: Estimates of tidal field strength, molecular gas densities and re- 
sulting predictions for t4, the tdis of a cluster with an initial Mc = lO'^ M a . Co l- 
umn 4: emperically derived valu es of tj are given, taken from: ^'Gicles et al. I ()2005l ): 

iLamers. Gieles &: Portegies ZwartI (|2005|): ^Lamers. Gieles. Bastian. et at. (,2005 '): *Parmentier 
fc de Griis (2007): iKrienke fc Hodgd (|2004f l: iBoutloukos fc LamersI (l2003h . 



Galaxy Tidal field Molecular gas density Predicted t4 Observed ti 

Rg/Vg [Myr] pn [10-^ Mq pc"^] [Gyr] [Gyr] 



M51^ 


10 


450 


0.13 


0.1 


M33^ 


15 


25 


1.4 


0.6 


Solar neighbourhood^ 


35 


30 


1.6 


1.3 


LMG* 


30 




< 6.6 


>1 


NGC6822'^ 


~ 35 




< 7.7 


rsj 4 


SMC^ 


40 


0.5 


8.2 


8 



The scaling of tjig with cluster density (Mc/r^) combined with the observed weak 
dependence of r^ on Mc, rh oc Mc°'^'^, results in a similar scaling of the mean t^^'^ with 



Mc as found for i^';-^ , i.e. oc Mc"'^ (fTTjl . 

For the solar neighbourhood (pn ~ 0.03 Mqpc"'^) t^;^^ = 2 Gyr, which combined 
with the tidal field (jl.ll) nicely explains the emperically derive d tdis of 1.3 Gyr and th e 
observed age distribution of clusters in the solar neighbourhood ( Lamers fc Giele^ 2006f ) . 

From (jl.ip and p.2p we see that the predicted idis scales with the tidal field strength 
(Rq/Vg) and the inverse of the molecular gas density (1/pn)- In table[T]we give values for 
these para meters for six galaxies, combined with predic t ions f or t^. The values for are 
taken from Gieles. Por tegies Zwart. Baumgardt. et all (l200 6l) (and references therein), 
Hever et al. ( 2004 ); iLeroy. Bolatto, Stanimirovi c et al\ (|2007) for the solar neighbour- 
hood, M51, M33 and the SMC, respectively. In the next section we compare this to 
empirically derived values of idis- 



2. Comparison to observations 

2.1. Empirically derived tdis values in different galaxies 
Under the assumption that tdis scales with Mc, iBoutloukos fc LamersI (|2003[ ) (BL03) 



introduced an empirical disruption law: tdis — ti (Mc/IO^M©)^. The value of ^4 and 7 
can be derived from the age and mass distributions (see BL03 for details). BL03 found a 
mean 7 of 7 = 0.62, agreeing nicely with ()1.1|) and (|1.2p . and values for ranging from 
^ 100 Myr to 8 Gyr. We summarise values of ^4 of clusters in six different galaxies 
taken from more recent literature in table [T] 

Note that Rq/Vq and l/pn roughly increase with increasing ^4. The variation in Rq/Vq 
is too small to explain the variation in ^4, which implies that in the galaxies with short ^4 
the disruption is dominated by GMC encounters. From Table\^we see that ttie decreasing 
trend in the emperical t/^ can be explained by increasing gas density and increasing tidal 
field strength. 

2.2. The clusters of the SMC 

A lot of attention has gone recently to the age distribution (dN/dt) of clusters in the SMC. 
Rafclski fc Zaritskv (2005) (RZ05) found that dN/dt is roughly declining as t~^, which 
Chandar et al.i (|2006( ) explain by mass independent cluster disruptiorU removing 90% 



f In fact the authors call their disruption model "infant mortality", but we prefer to reserve 
this term for the dissolution of clusters due to gas expulsion. In addition, 3 Gyr old clusters 
have survived 25% of a Hubble time, so they are not really infant anymore. 
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Figure 2. The age distribution (dN/dt) of star clus t ers in the SMC as found in different studies 
in literature. The data se t of [Rafelski fc Zaritskvl ll2005f) (R Z05) is very incomplete for low 
mass clusters at old ages l|Gieles. Lamers fc Portegies Zwart| [2007'). so a mass cut at 10^ Mq 
was applied. The general trend found in these studies is that dN/dt is flat up to an age of 
0.3 — 1 X 10^ yr and then it declines as t~^''' . The dashed line is the predicted slope for dN/dt 
at old ages when tdis oc Mc^, with 7 = 0.62. 



of the clusters each age dex. iGieles. Lamers fc Portegies ZwartI (j2007l ) showed that the 
decline is caused by incompleteness and that the AN/dt is flat in the first ^ 1 Gyr when 
using a mass limited sample. The dN/dt based on ages which are derived from extinction 
corrected colours starts declining a bit earlier than the one based on uncorrected colours 
(figure [2]). However, the general shape is similar to tha t found by other authors: a flat 
part in the flrst 0.3 — 1.0 Gyr (recently reconfirmed bv Ide Grijs fc GoodwinI 2007 ) and 
then a steep decline (oc t~^-''). When tdis oc M^'^ , then the dN/dt at old ages declines as 
i-i/^ for both mass and magnitude limited samples (BLG3). The decline of t implies 
7 ~ 0.6, in agreement with the theoretical predictions p. II and ll.2p . 



2.3. Selection effects and biases: a cautionary note 

Observed cluster samples are always heavily affected by the detection limit, causing the 
minimum observable cluster mass (Mmin) to increase with age, due to the fading of 
clusters. To illustrate this effect we create an artificial cluster population with a constant 
cluster formation rate (CFR) and with a power-law CIMF with index —2. In the left 
panels of figure [3] we show the ages and masses (bottom) and the corresponding dN / dt 
(top) when the sample is mass limited. The dN/dt is flat which is the result of the 
constant CFR we put it. In the right panel we remove the clusters which are fainter than 
Mv — —4.5. The mass of a cluster at the detection limit, Afmin(i), increases with age 
as Afmin(i) oc 0.4M^SP(t), where M'^^^{t) is the evolution of My with age from an SSP 
model. For a power-law CIMF with index —2, the resulting dN/dt scales with Mmin as 
dN/dt (X l/A/min(t) (BL03), which is shown in the top right panel of figure [3] 

The detection limit is usually expressed in My. However, deriving cluster ages from 
broad band photometry requires the presence of blue filters such as U and B. We show 
the M,nin(i) for a [/-band detection limit of Mu = — 5 as a dashed line in the age vs. 
mass diagram. The resulting dN/dt (shown as a dashed line in the top right panel) 
declines approximately as dN/dt oc t~'^'^, i.e. steeper than the 1/^-band prediction. It is 
of vital importance to understand the effect of incompleteness in different filters before a 
disruption analyses can be done based on the slope of the dN/ dt distribution. 
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Simulated cluster population: constant CFR, no disruption 
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Figure 3. Simulated ages and masses of a cluster population that has formed with a constant 
cluster formation rate (CFR) and with a power-law CIMF (A*' oc M~^). In the left panels we 
show the result of mass limited sample, with Afjim — 10 Mq . In the right panels we assume that 
the sample is magnitude limited, with Mv,Uin = —4.5. The limiting mass due to a magnitude 
limit and the resulting prediction for dN/dt of a magnitude limited sample are shown as full 
lines (red). The prediction for a [/-band limit {Mu = — 5) is shown as dashed lines (blue). The 
dotted lines show power-law approximations for the predicted shapes of dN/dt. 
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